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We estimate the Collins fragmentation function by introducing the effect of gluon rescattering in 
a model calculation of the fragmentation process. We include all necessary diagrams to the one-loop 
level and compute the nontrivial phases giving rise to the Collins function. We compare our results 
to the ones obtained from pion rescattering. We conclude that three out of four one-loop diagrams 
give sizeable contributions to the Collins function, and that the effect of gluon rescattering has a 
magnitude comparable to that of pion rescattering, but has opposite sign. 
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I. INTRODUCTION 

The Collins fragmentation function [jj] is an example of how the orientation of the quark spin can influence the 
direction of emission of hadrons in the fragmentation process. The existence and the features of this function are 
related to important questions of nonperturbative QCD, such as the role of chiral symmetry and color gauge invariance 
O • in the hadronization process. The Collins function is believed to be at the origin of single-spin asymmetries in hard 
hadronic reactions 0, IE 0, IE IE > which lead to attempts of estimating it from phenomenology [E IE [HI EH d] • 
■ The Collins function is a so-called T-odd entity. T-odd functions typically require the interference between two am- 
' plitudes with different imaginary parts to exist. Perturbative calculations of fragmentation functions in quantum field 
theories can provide - through loop corrections - the necessary nontrivial phases in the fragmentation amplitude . 

The Collins fragmentation function has been estimated in a chiral invariant approach where the effective degrees 
of freedom are constituent quarks and pions, coupled via a pseudovector interaction In order to generate the 
required phases, one-loop corrections to the quark propagator and vertex have been included. 

In the meanwhile, gluon loop corrections in distribution functions have been investigated. In the context of a 
spectator model of the nucleon, it has been shown that the exch ang e of a gluon between the struck quark and the 
target spectators gives rise to T-odd distribution functions | l5l llq . This has been interpreted as the one-gluon 
£T) . approximation to the gauge li nk, which is included in the definition of parton distribution functions and insures their 
color gauge invariance flfi lr7L llSl . 1131 . It was soon realized that gauge link effects have different signs in semi-inclusive 
DIS and Drell-Yan scattering 0,HE|, implying different signs for T-odd distribution functions in the two processes. 
This is for the first time an example of violation of universality of the distribution functions, though it consists simply 
of a sign change. Gluon rescattering has been used thereafter to estimate T-odd distribution functions |2lLl2^|. 

As in the distribution functions, the gauge link can generate nontrivial phases and T-odd effects in the fragmentation 
functions as well. This was analyzed in Ref. pflj . where it was shown in particular that one-gluon contributions to 
• • . the gauge link don't change sign in semi-inclusive DIS and in e + e~ annihilation. 

In this work, we model the tree-level fragmentation process in the same way as in Ref. ,14'j . To generate T-odd effects, 
instead of introducing pion-loop corrections, we include gluon rescattering corrections. We take into consideration 
all possible one-loop amplitudes, including those that arise from the gauge link. We compute the Collins function 
and its first two moments, which can be experimentally accessed in single-spin asymmetries in semi- inclusive DIS, as 
well as in e + e~ annihilation. A similar calculation has been presented very recently by Gamberg et al. [23|, using a 
pseudoscalar coupling complemented with a Gaussian form factor to model the tree level fragmentation. However, in 
that work only one of the four possible gluon rescattering contributions has been taken into account. 
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FIG. 1: Single gluon-loop corrections to the fragmentation of a quark into a pion contributing to the Collins function. 
Hermitian conjugate diagrams (H. c.) are not shown explicitly. 
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II. CALCULATION OF THE COLLINS FUNCTION 



We use the following definition of the Collins function [2^, H 

F ij k 



Ht{z, z*K) = Tr [A(«, K) 1^-75] . 



(1) 



The correlation function A(z, k T ) can be written as [l? 



A{z,k T ) 
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4z(2tt) 3 
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where the notation £r w indicates a straight gauge link running from a to b. In this work, we will use the Feynman 
gauge. In the case of transverse-momentum independent fragmentation functions, i.e., after integrating over k T , by 
choosing a light-cone gauge the link can be reduced to unity. However, in the case of transverse-momentum dependent 
functions - as the Collins function - the gauge link cannot be neglected and becomes one of the possible sources of 
nontrivial phases in the fragmentation amplitude and thus of T-odd fragmentation functions . Note that in case 
of transverse-momentum dependent functions with lightlike Wilson lines, divergencies can appear |27l l28j . However, 
for our study this problem does not show up. 

To model the tree-level fragmentation process, we make use of the chiral invariant effective model of Manohar and 
Georgi [29| . The unpolarized fragmentation function in this model reads 0] 



D 1 (z, 
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To obtain a nonzero Collins function, we have to compute one-loop contributions. In contrast to what was done 
in our previous calculation |l4j . instead of introducing pion loops, we now consider the effect of including gluon 
rescattering. 

In Fig. ^ we illustrate the Feynman diagrams involved in the calculation of the Collins function. Diagrams (a) and 
(b) - the self-energy and pion vertex corrections - have analogous contributions in the pion-loop case. However, due 
to gauge invariance we have to take into consideration also diagrams (c) and (d), which we will henceforth call the box 
and photon vertex corrections. Despite their appearance, these contributions do not break factorization, since they 
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FIG. 2: One-loop self-energy, pion vertex, box, and photon vertex corrections. 



can be interpreted as contributions to the gauge link, which is included in Eq. J2J lfiL UtI ITsl ITflj . In fact, diagram 
(d) contributes not only to the gauge link of the fragmentation function. However, since in this specific calculation 
its contribution turns out to vanish, we will avoid considering this issue. 

In all these diagrams, the gluon loops can give rise to nontrivial phases to the scattering amplitude. Through the 
interference with the tree-level amplitude, these contributions generate nonzero T-odd fragmentation functions, such 
as the Collins function. The portions of the diagrams relevant for the generation of imaginary parts are sketched in 
Fig. [2] They can be expressed analytically as 
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where Cf — 4/3. We choose the kinematics in the following way (in light-cone coordinates): 
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To identify the contributions to the Collins function, we write down the explicit expressions for the cut diagrams 
of Fig. n and we apply the definition of the Collins function given in Eq. (JJJ. 

It turns out that only some specific elements of the imaginary parts of the diagrams in Fig. contribute to the 
Collins function. For simplicity we will denote them as Imcr, Ini7, Im£, Im^. The formula of the Collins function 
can be written then as 
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Considering only the leading power in Q and assuming the gluons to be massless, the calculation of the four diagrams 
yields 
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where the integrals introduced above correspond to 
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In the previous formulae, we made use of the Kallen function, X(k 2 ,m 2 ,m 2 ) — [k 2 — (to + rn 7r ) 2 ][fc 2 — (m — m,) 5 
Note that the final result for the Collins function is independent of Q 2 . 



III. NUMERICAL RESULTS 



We perform the numerical integration over the transverse momentum according to 



Di(z)=it 



dK 2 Di(z,K£), 



(18) 



o 



where K T = —zk T . We impose a cutoff on the virtuality of the fragmenting quark, so that k 2 < /i 2 . Due to the 
kinematics, this choice fixes the upper limit of the K 2 integration to 



^Tmax = z C 1 - z ) V 2 - z m2 - (1 ~ z ) ml . 



For the numerical computations, we use the following values of the parameters of our model: 



m q = 0.3 GeV, g A = 1, /i 2 = 1 GeV 2 



a, = 0.3. 



(19) 



(20) 



The choice of the first three parameters has been extensively discussed in Ref. . The choice of the strong coupling 
constant a s corresponds to a reasonable value for Q 2 ss 1 GeV 2 . 
In Fig. |3 and 0| we plot the ratio 
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FIG. 3: Result for _ff 1 L ' 1,/2 - ) / D\ including only the self-energy and vertex gluon-loop corrections (solid line) and comparison 
with the result of the pion-loop model (dashed line) from Ref. [l4j . 
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FIG. 4: Result for jD\ including all gluon-loop contributions (solid line) and comparison with the result of the pion-loop 

model (dashed line) from Ref. [14j. 
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FIG. 5: Result for H^ 1 ^ /D\ including all gluon-loop contributions (solid line) and comparison with the result of the pion-loop 
model (dashed line) from Ref. [l4]|. 



which enters the unweighted transverse single spin asymmetries for pion production in semi-inclusive DIS BUI. In 
Fig. [3|we show only the contributions of the self-energy and pion vertex corrections, diagrams (a) and (b) of Fig.^ 
They have a direct correspondence to the pion-loop case. As shown by the plot, these contributions are smaller and 
have an opposite sign compared to the pion-loop ones. In Fig. we show the sum of all gluon-loop contributions, 
i.e. all diagrams of Fig. ^ including in particular the box diagram, which contributes to the gauge link and has no 
analogous term in the pion case. 
In Fig. we plot the ratio 



H 



-L(l) 



Ht(z,K 2 T ), 



Di(z) D x (z) J " T 2z 2 ml 
which typically enters weighted transverse spin asymmetries in semi- inclusive DIS |l4l l30t l3l| . 



(22) 



IV. CONCLUSIONS 



In this work we have calculated the Collins function for pions describing the fragmentation process at tree level 
by means of a chiral invariant effective theory, and generating the required nontrivial phases by means of gluon 
rescattering. We have computed all necessary diagrams contributing to the Collins function at the one-loop level, 
including the ones involved in the gauge link. Out of the four diagrams we considered, one gives no contribution, 
while the other three have similar magnitudes and cannot be neglected. A word of caution is in order at this point: 
we have performed a perturbative expansion in as in a regime where its justification can be questioned. Moreover, 
we emphasize once again that this is only a model calculation, whose reliability might be very limited. Nevertheless, 
we believe that it could give an indication of the influence of the gluon dynamics on the Collins function. 

We compared our results with those formerly obtained in Ref. 14), where pion rescattering was considered as a 
possible source of nontrivial phases. The two approaches are identical at tree level, allowing a clear comparison of the 



differences at the one-loop level. We have found that the gluon rescattering mechanism produces a Collins function 
with opposite sign compared to the pion rescattering mechanism. The two effects are similar in magnitude, except 
perhaps at high z, where pion rescattering dominates. Cancellations between these two competing mechanism could 
decrease the experimental asymmetries. 
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